Gradient estimate of an eigenfunction on a 
compact Riemannian manifold without boundary 
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Abstract. Let e^(x) be an eigenfunction with respect to the Laplace-Beltrami operator Am on a 
compact Riemannian manifold M without boundary: Am£\ = A 2 e^. We show the following gradient 
estimate of e\: for every A > 1 , there holds A|| eA || oo /C < || VeA ||oo < CA|| ||oo , where C is a positive 
constant depending only on M. 
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1 Introduction 

> 

\Q ' Let (M, g) be an n-dimensional compact smooth Riemannian manifold without boundary 

and Am the positive Laplace-Beltrami operator on M. Let L Z (M) be the space of square 
integrable functions on M with respect to the Riemannian density dv(M) = \/g(x) dx := 
^/det~{g^) dx. Let z-\ (x), ezi x )> " " " be a complete orthonormal basis in L 2 (M) for the 
eigenfunctions of Am such that = Ao < A 2 < A2 < • • • for the corresponding eigenvalues, 
where ej(x) (j = 1,2,---) are real valued smooth function on M and Aj are nonnegative real 
numers. Also, let ej denote the projection onto the 1-dimensional space Ce^. Thus , an L 2 
function f can be written as f = ^^ ej (f ), where the partial sum converges in the L 2 norm. 
Let A be a positive real number > 1. We define the spectral function e(x,y,A) and the unit 
band spectral projection operator xa as follows: 

e(x,y,A) := Y_ ^jM^iv) , 

A, <A 
Aj£(A,A+l] 

In 1968, Hormander [6] obtained a one-term asymptotic expansion of the spectral func- 
tion of a positive definite elliptic linear operator, whose Laplacian case, also called by the 
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local Weyl law, has the expression as follows: 

e(x,x,A) =A n /((47t) n/2 r(l+n/2))+0(A n - 1 ), A^oo. 
As a consequence Hormander proved the uniform estimate of eigenfunctions for all x G M: 

Y_ |e j (x)| 2 <CA- 1 . 

Aj e(A, A+1] 

We note here that in the whole of this paper C denotes a positive constant which depends 
only on M and may take different values at different places. C. D. Sogge noted in (1.7) of 
[12] the equivalence of the above estimate with the following L°° estimate of xa, 

IIXAflloo < CA (n -^ /2 ||f|| 2 . (1) 

Here |jf|| r (2 < r < oo) means the L T norm of the function f on M. As noted in Lemma 2.7 
by the second author [14], there also holds the similar equivalent relationship between the 
following two gradient estimates for eigenfunctions and xa: 

Y_ |Vej(x)| 2 < CA n+1 forallxeM, (2) 

A, g(A,A+1] 

and 

IIVxAflloo < CA<- +1 »/ 2 ||f|| 2 . (3) 

Here V is the Levi-Civita connection on M. In particular, Vf = Y_j g l '9f/9xj is the gradient 
vector field of a C 1 function f, the square of whose length equals Y.\ j 9 v '(3"f/3xi)(3"f/3x-j), 
where (g 1 ') is the inverse of the metric matrix gy. By using the wave group, Yu Safarov 
and D. Vassiliev proved a very general theorem (Theorems 1.8.5 and 1.8.5 in [9]) on the 
spectral function of a positive definite elliptic linear operator so that the gradient estimate 
(2) is its immediate corollary. By using the parametix of the wave operator 3 2 /3t 2 + Ajw, 
the second author [14] also proved in a slightly different way the Laplacian case of the afore- 
mentioned theorems of Yu Safarov and D. Vassiliev, which is also sufficient to deduce the 
estimate (2). X. Xu [17] applied the maximum principle argument to proving (3), and used 
this estimate to give a new proof for the Hormander multiplier theorem. A. Seeger and C. 
D. Sogge [11] firstly proved the the Hormander multiplier theorem by using the parametrix 
of the wave kernel. In this paper we stick at estimating the gradient of a single eigenfunction. 

Theorem Let e\{x) be an eigenfunction with respect to the positive Laplace- Beltrami 
operator Am on an n- dimensional compact smooth Riemannian manifold M without 
boundary: AM.^\ = '\ 2 e\. Then, for every A> 1, there holds 

M\ e \\\oo/C < HVexHoo < CA||e A ||oo, (4) 

where C is a positive constant depending only on M. 
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We can also obtain the above estimates (2) and (3) by Hormander's local Weyl law and 
the same argument of Theorem. The details will be given in Section 4. The second author 
announced a more general estimate for ||V k e A ||oo with k a positve integers in Theorem 5.1 
of [15]. However, the method we use in this paper can not be applied to deducing higher 
derivative estimates. We plan to discuss this question in a future paper. 

We conclude the introduction by explaining the organization of this paper. Sections 2-4 
contain the proof of Theorem. In Section 2, we show the lower bound of the gradient by 
the equidistributional property of the nodal set of eigenfunctions. In Section 3 we reduce 
by rescaling the gradient estimate from above to that of an almost harmonic function. In 
Section 4, we prove the gradient estimate from above for the almost harmonic function by 
Yau's gradient estimate and the Green function for compact manifolds with boundary. We 
also in this section sketch a proof of the above estimates (2) and (3) by the same argument 
with Theorem, which is different from those in [17] and [14]. In the last section we propose 
a conjecture that a similar gradient estimate as Theorem holds for a compact Riemannian 
manifold with boundary. 



2 Nodal sets of eigenfunctions 

The nodal set of an eigenfunction e A is the zero set 

Z £A :={x£ M:e A (x) =0}. 

A connected component of the open set M\Z eA is called a nodal domain of the eigenfunc- 
tion e A . The following fact, due to J. Bruning [2], is critical in deducing the lower bound 
estimate. The reader can also find a proof in English in Theorem 4.1 of S. Zelditch [19]. 

Fact 1 Under the assumption of Theorem, there exists a constant C such that each 
geodesic ball of radius C/A in M must intersect the nodal set Z eA of e A . 

Proof of the lower bound part of Theorem Take a point x in M such that |e A (x)| = 
1 1 e A | |oo- Then by Fact 1, there exists a point y in the ball B(x, C/A) with center x and radius 
C/A such that e A (y) =0. We may assume A so large that there exists a geodesic normal 
chart (r,9) G [0, C/A] x § n_1 (1) in the ball B(x, C/A). By the mean value theorem, there 
exists a point z on the geodesic segment connecting x and y such that 

A A 



3r 



3 



3 Eigenfunctions on the wavelength scale 

In this section we review quickly the following principle: 

On a small scale comparable to the wavelength 1 /A, the eigenfunction e A behaves like 
a harmonic function. 



This principle was developed in H. Donnelly and C. Fefferman [3] [4] and N. S. Nadirashvili 
[8] and was used extensively there. Our setting of the principle in the following is borrowed 
from D. Mangoubi [7], where he applied it to studying the geometry of nodal domains of 
eigenfunctions. 

Fix an arbitrary point p in M and choose a geodesic normal chart (B(p,6), x = (xi,- ■■ ,x n )) 
with center p and radius 6 > depending only on M. In this chart, we may identify the ball 
B(p,6) with the n-dimensional Euclidean ball B(6) centered at the origin 0, and think of the 
eigenfunction e A as a function in B(6). Our aim is to show the inequality |Ve A (p)| < CAHe^. 
Actually we plan to prove a slightly stronger one 

3e Ai 



OXj 



< CA||e A (x)|| L oo( B ( 1/A))) (5) 



which we remember that g-ij(O) = 6^ in the normal chart. The rough idea is to consider 
the rescaled function cb A (y) = e^Oy/A) in the ball B(1) instead of the restriction of the 
eigenfunction e A (x) to the ball B(1/A), where we assume A so large that 1/A < 6/2. The 
above estimate is equivalent to its rescaled version 

9cb A 



, (0) 
3yj 



< C||cp A (y)|| L oo (B(i))- (6) 

. " »N i 

On the other hand, rescaling the eigenfunction equation Am^a = ^ 2e A tells us that the 
function cb A behaves like a harmonic function. The details are given in what follows. 
The eigenfunction equation Ajvie A = A 2 e A in B(1/A) can be written as 

1 



— 7= V 3 Xl (g v ^d Xj e x ) =A 2 e A . 
V9 y 



Hence the function cb A satisfies the rescaled equation in B(1) 

1 



ll3 m (g^9 yj cb A )=cb A , (7) 

where gij )A (y) = gij(y/A), g%[y) = g 15 (y/A) and ^/gT(y) = (^g)(y/A). We note that the 
above equation coincides with 

A A cb A = cb A) 

where A A is the positive Laplace-Beltrami operator with respect to the metric gy iA on 
B(l). We observe that the scaling xhij = Ax gives an isometry from (B(1/A), A 2 gtj) onto 
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(B(1), gij ; a) - Therefore, the open Riemannian manifold (B(1), gij t \) has uniformly bounded 
sectional curvature for all A > 1 . This fact will be crucial in next section. 



4 Estimates from above 

We prove (6) in this section by using the notations in Section 3. 
Step 1 Recall that Aac|)a = <t>A in B(1 ). We can write the function 4>a as the sum §\ = 
ua + va such that the functions ua and va satisfy the following boundary- value problems: 

f A A v A = 4>A in 1(1), 
\ v A = on3B(l), 

Step 2 By the gradient estimate for harmonic functions in page 21 of R. Schoen and 
S.-T. Yau [10] and the maximum principle, there holds 




3yj 



< C||uA(y)||i°°(B(i)) < C||cpA(y ) ||t_°° (B(i))- 



Here we also use the fact that the Ricci curvature tensor of the metrics g^A are uniformly 
bounded for all A > 1 . 



Step 3 Let G(y,z) be the Green function on the compact manifold (B(1), gij,\) with 
smooth boundary such that in the sense of distribution AA, z G(y, z) equals the Dirac delta 
function 6 y (z) at y and G(y,z) vanishes whenver y or z belongs to the boundary 3B(1). 
Remember that the sectional curvatures of gy^ are uniformly bounded for all A > 1 . Care- 
fully checking the explicit construction of the Green kernel in pages 106-113 in T. Aubin 
[1], we find that 

KVyGKCz^CdtCz) 1 -- 

Taking gradient at in the Green formula for va 



G(y,z)cMz) dV(z) 



gives the estimate |(Vva)(0)| < CIIc^aIIl 00 (B(1))- Remember that gy^lO) = 6y so that |(Vva)(0)| 
is comparable to Y.) l(3 v A/3yj)(0)l- The proof is completed. 



Following the above arguments, we conclude this section by sketching an alternative 
proof of the estimates (2) and (3), which were proved by X. Xu [17] and the second author 
[14] by using the maximum principle and the parametrix of the wave kernel, respectively. 
We also use the notations above. Take a square intergable function f on M and consider the 
smooth function u:=xx^, and the Poisson equation v := Ajviu in the geodesic normal chart 
B(l/A) centered at a given point p in M and with radius 1/A. Rescaling the ball B(l/A) 
with metric gy to B(l) with metric g^A, and using the L°° estimate 

||Xxf||i~ < CA^-^Hflb 
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implied by the local Weyl law and the Cauchy-Schwarz inequality, we reduce the question 
to showing 

< C||u|| L oo (B(1/A)) + CA 2 ||va||l oo (B(1))> 

where u A (y) = u(y/A) and v A (y) = v(y/A). But the last inequality follows from the argu- 
ments in Steps 1-3 above. The proof is completed. 

5 Compact manifold with boundary 

D. Grieser [5] and C. D. Sogge [13] obtained a similar sup norm estimate as (1) for xa 
associated with either the Dirichlet or the Neumann Laplacian on a compact Riemannian 
manifold N with boundary. X. Xu [16] [18] obtained similar gradient estimates as (3) for 
both the Dirichlet and the Neumann boundary value problems on N by using a clever 
maximum principle and the results of Grieser and Sogge. However, we could not use the 
maximum principle argument in their papers to deduce more refined gradient estimate for 
a single eigenfunction on N. We think that new ideas should be introduced to resolve the 
following conjecture. 

Conjecture Let N be a compact Riemannian manifold with smooth boundary and 
e A (x) be an eigenfunction of either the Dirichlet or the Neumann Laplacian An on N: 
Ai\je A = A 2 e A . Then there exists a positive constant C depending only on N such that 
for all A > 1 there holds 

A||e A ||oo/C < || V e A ||oo < CA||e A ||oo. 
We should point out that the lower bound estimate for the Dirichlet eigenfunction ex, 

||Ve A ||oo > CA||e A || 

follows easily from the similar argument in Section 2. We also observe that the upper bound 
estimate |Ve A (x)| < CAHe^Hoo for every x outside the boundary layer {x G N : d(x, 3N) < C/A} 
follows from the argument in Sections 3 and 4. The authors will discuss this question in a 
future paper. 
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